Depinning transitions in elastic strings 
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We study the depinning transitions of elastic strings in disordered media in two different cases. We 
consider the elastic forces to be of infinite range in one case, where the magnitude is proportional 
to the extension of the string. The critical points and elastic behavior can be estimated to some 
extent in this case. The exponents are estimated numerically and scaling relations are found to be 
obeyed. We have also considered a model where the elastic force is constant in magnitude. The 
critical points can again be argued and the critical exponents are numerically estimated. The scaling 
relations are well obeyed. Both of these models do not fall into known universality classes. 



I. 



INTRODUCTION 



Dynamics of elastic strings, driven through a medium 
with quenched disorder has served for modelling of 
the dynamics of diverse physical phenomena such as 
charge density waves 0, vortices in type-II supercon- 
ductors Q, magnetic domain walls in external magnetic 
field [i[ and also propagating fracture fronts 

H (for a 

recent review see [6|). 

When an elastic string is driven through a random 
medium, the forces that act on the string are the follow- 
ing: the uniform external force (/), the random pinning 
force due to impurities and the local or non-local force 
due to elastic stretching of the string (f e i)- When the 
external driving force is small enough, the string as a 
whole would not propagate (in the long time limit). In 
that case the system is said to be pinned by the impu- 
rities present in the medium. On the other hand, when 
the force exceeds a threshold f c , the string will propagate 
through the disordered medium with a finite velocity (in 
the long time). At / = f c , string is just depinned. 

Global failure processes due to fracture, has been stud- 
ied using the tools of statistical mechanics before Q. 
However, a long studied approach to study depinning 
transition problems have been to deal the depinning of 
the string as a dynamical phase transition. Much atten- 
tion have been paid in this area both theoretically as well 
as experimentally (see e.g., [8rll3j). Since in the depinned 
state the average velocity v of the string is non-zero and 
in the pinned phase it is zero, the velocity is taken as the 
order parameter for this transition. As it is a second or- 
der phase transition, the order parameter should scale as 
v ~ (.f ~ fc) e where is velocity or order parameter expo- 
nent. Other critical exponents and scaling relations can 
also be derived to formally classify a universality class of 
this transition. 

The most studied universality classes in the depinning 
transitions are the Edwards- Wilkinson (EW) |14| - fl6| and 
Kardar-Parisi-Zhang (KPZ) [l7l - [l9} universality classes 



(see Ref.[l[ for details). Apart from these models where 
the elastic force is essentially local (depends on the po- 
sitions of the nearest neighbors of an element), there are 
some non-local extensions, which are relevant for fracture 
front propagation, where the elastic force on an element 
has contributions from all other elements. However, the 
amount of contribution decreases in a power law (mostly 
inverse square) with the distance from the said element 
(see e.g., |2(j). A mean field model, in which the elas- 
tic force on an element depends on its distance from the 
average height, has also been studied (2lj . 

In this paper we study two different limits of the form 
of the elastic force. First we study the case where the 
elastic force is proportional to the total elongation of 
the string and has equal magnitude at each point on the 
string. The sign, however, depends upon the sign of the 
local curvature. Clearly, in this case, the elastic force is 
of infinite range. Using this fact, we argue to find the 
critical point at which the string is just depinned. The 
critical behavior is obtained by integrating the equation 
of motion. The exponent values satisfy usual scaling re- 
lations. 

Secondly, we consider where the magnitude of 

the elastic force is a constant and the sign depends upon 
the local curvature. Furthermore, the pinning forces re- 
main constant for a given point on the string. This is 
in contrast with the convention that the pinning force 
changes its value after certain intervals along the propa- 
gation path of the string. 

The paper is organised as follows: In the next section 
we describe the model and the numerical methods used 
in obtaining the height profile. Next we present the scal- 
ing analysis and estimate the critical exponents and the 
scaling laws satisfied by them. Finally we discuss these 
results and conclude. 



II. 



MODELS AND METHODS 
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As mentioned before, an elastic chain driven through 
a medium with quenched disorder faces three types of 
forces: External driving force (/), elastic force (f e i), 
which tends to conserve the length of the string and 
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a random pinning force (f p ) due to quenched disorder. 
The (discretized) equation for the time evolution of the 
"height" profile of the string is 



hi{t + l) = hi{t)+f el + f p + f 



(1) 



where, hi(t) denotes the height of the i-th element of the 
string at time t. The movement of the string is allowed 
only in the forward direction. 




FIG. 1: The height profile hi(t) is plotted against i for 

different time instants for EW, KPZ and the two 
models studied here. The time index increases upwards. 

One can define an average velocity (at time t) of the 
string as 



v(t) 



1 N 



hi(t-l)). 



(2) 



i=i 



The 'pinned' and 'depinned' phases of the string is de- 
termined from the fact whether the above quantity van- 
ishes or attains a finite value in the long time limit. 
This is analogous to non-equilibrium active-absorbing 
type phase transitions with velocity (v(t)) playing the 
role of the order parameter. The external force / is usu- 
ally taken as the 'ordering field'. Therefore, near the 
transition point one would expect 



v(t oo) ~ (/ - f c y 



(3) 



where, f c is the critical value of the external force where 
the chain is just depinned. Other critical exponents and 
scaling relations (discussed later) can also be obtained, 
which define the universality class of the transition. The 
universality classes are dependent upon what form of the 
clastic force is chosen. For example, when the elastic 
force is simply proportional to the local curvature (V 2 /i), 
the universality class is Edwards- Wilkinson (EW); when 
it is proportional to the extension of the chain in the 
nearest neighbor, it is KPZ type and so on. All these 
cases are for uncorrelated noise, the cases of correlated 
noise change the situation a lot (see [l|) which we shall 
not discuss here. 



Below we study two models, where this elastic forces 
are different. In one model, the elastic force is infinite 
range and its magnitude is proportional to the total ex- 
tension of the string and the sign of the force depends 
upon the local curvature. In another model, the mag- 
nitude is taken as a constant and the sign still depends 
on the sign of the local curvature (Fig. [TJ depicts height 
profiles of these two and EW and KPZ models at critical 
point for different time instances). 



A. Model I 

Here we consider an elastic chain of infinite range inter- 
action driven through a medium with quenched disorder. 
We assume linear elastic behavior, such that the net ex- 
tension of the chain generates an elastic force which is felt 
by all element of the chain. In particular, we consider the 
following equation of motion for the height profile 



hi(t + 1) = hi(t) + f e i + rn(hi(t)) + f 
L ^ 



hi(t) 



(M*)-^+iW) 2 + i-i 



+m(hi(t)) + f 



Ci 
(4) 



^/(h 3 (t)-h J+1 (t)y 



/L, Ci 



with \f e i\ = 

j 

sgn(hi+i + hi-i — 2hi) and hi(t) denotes the height of the 
i-th element of the discretized chain at time t, rji denotes 
the random uncorrelated quenched noise in the pinning 
force (the nature of which we further illustrate below), / 
denotes the constant external force applied uniformly (on 
each element of the chain) and Ci denotes the sign of the 
local curvature i.e., the net extension force acts to restore 
the flat initial condition for the chain thereby creating 
a restoring force for the locally advanced elements and 
pushes the locally behind elements. 

The nature of the disorder in pinning force rji is such 
that with probability p it is zero and with probability 
1 — p it is uniform between and —a. p acts here like a 
probability of presence of vacancies within the material 
(in the following we will study p values ranging from 
to 0.5). To take the pinning strength of the material to 
have uniform values is a matter of convention, the critical 
behavior does not change if one takes it differently (trian- 
gular, say). It is to be noted here is that, unlike general 
convention, the disorder rji here has no positive value. 
In our simulations, we took an elastic chain of length L 
{— 5000 mostly), with periodic boundary conditions. 



1. Integration of the equation of motion 

The critical properties (exponent values) of this transi- 
tion are to be found out by numerically integrating the 
equation of motion (see e.g., [H[). This is done following 
the equation: 
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hi(t + At)= hi(t) + Gi(t)At if Gi(t) > 

= hi(t) otherwise (5) 

where Gi{t) = f e iCi + r]i(hi(t)) + f and hi{t) denotes 
the integral part of hi(t), At is a small interval of time, 
mostly taken as 0.01 but our results do not change even 
if it was made smaller (say, 0.005). Clearly, it is ensured 
in the dynamics that any element never takes backward 
step, which is unphysical from the point of view of frac- 
ture front propagation. 

In the following, we have presented numerical results 
for L — 5000, St = 0.01, p = 0.5 and averaging over 100 
configurations. It is checked that making St = 0.005 or 
L = 10 4 does not change the conclusions. 

B. Model II 

In this case we make two major changes. First, the elas- 
tic force is taken as a constant and its magnitude is de- 
termined by the local curvature. In this sense, it is a 
model with short range interaction. Secondly, we keep 
the quenched disorder for a given site fixed in height, 
i.e., rji is no longer a function of (integral part of) hi(t). 
The evolution equation for height, therefore, becomes 

h i (t + l) = h i (t)+FC + r ]i + f (6) 

where F is the constant magnitude of the elastic force 
and other symbols have same meanings as before. 

One can apply this model for 'depinning' transition in 
a flock of birds flying in a line (see e.g., [24| for refer- 
ences related to bird flocking). The elastic force denotes 
the tendency of a bird to fly along with its two nearest 
neighbors, the pinning force may be considered as some 
measure of '(un)fitness' of a bird (that is why it is taken 
as independent of height or in this case, distance trav- 
elled), and the external force is the urge or necessity of 
flying (lack of food, presence of enemy etc.). Of course, 
v — } limit does not apply for a flock of flying birds. But 
here we can measure the velocity from the rest frame of 
the slowest moving bird. Then the 'depinning' transition 
would indicate whether the birds fly 'together' or gets 
scattered in the long time limit. If the average velocity 
of the flock is zero with respect to the slowest moving 
bird, then they fly 'together' or the flock is 'pinned'. If, 
however, the average velocity is finite in the rest frame 
of the slowest bird, then in the long time the flock will 
be 'scattered' or the flock is said to be 'depinncd'. 

The nature of disorder used has similar properties as 
before. It has a random value between [0 : —1] with 
probability 1 — p and zero otherwise. Below we present 
results for p = case, but the exponent values do not 
change for finite p values. 

The evolution rule (Eq.([5])) is analogous to the equa- 
tion of motion presented above for an elastic chain. But 



the major difference is that in this case it is not an 'equa- 
tion of motion' but only a rule chosen for the movement 
of the birds. Therefore the choice of At is arbitrary and 
we fix it to unity in this case. 

C. Scaling analysis 

Using the numerical method described above we measure 
the height profile of the chain with time. At and near 
the depinning transition point (/ = / c ), several quan- 
tities related to the height profile show interesting be- 
havior. To begin with, we measure the average of the 
velocities of the elements, which is the order parameter 
v(t) = d{hi(t))i/dt. Above the depinning point, the av- 
erage velocity initially decreases, but at long time satu- 
rates to a finite non-zero value depending on the external 
force. On the other hand if the force is below the depin- 
ning point, the average velocity decreases faster than a 
power law and eventually becomes zero. Just at the crit- 
ical point the velocity decreases as a power law 

v(t) ~ r s . (7) 

Not only does this study give the value of the exponent 
S, it is also useful in locating the critical point. 

Considering the long term saturation values above the 
depinning transition, one would expect the growth 

v(t^^J)~(f-f c f : (8) 

where 9 is the velocity exponent. 

Another interesting quantity is the mean square fluc- 
tuation or the width of the surface, which is defined as 

W(L,t) = (i J2 Mt) -Ht))f) 1/2 - (9) 

i 

This will show a growth of the form W(t) ~ t" . 

Also it is argued that at criticality, the velocity of 
the chain is essentially governed by avalanches and an 
avalanche takes a time t to advance a distance equal to its 
width [2(j. Therefore, one would expect v(t) ~ W(t)/t, 
which would imply S + f3 = 1. 

Finally, the off-critical scaling of the average velocity 
has the form 

v(f,t)~r 5 F(t\f-f c n, (io) 

where v is the exponent for the diverging correlation 
time. It is expected to satisfy v = 9/5. 

Below we present numerical results for the two models 
estimating the above mentioned critical exponent values 
and show that the scaling relations are obeyed. 

III. RESULTS 

The critical points and for that matter the phase dia- 
grams for both the models can be argued to some extent. 
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However, the critical behavior, i.e., the values of the ex- 
ponents are to be found out numerically as indicated be- 
fore. In this section we present the phase diagrams and 
critical behaviors for both of the models. 



A. Model I 

1. Elastic behavior and critical points 

Due to the infinite nature of the elastic force, it is possible 
to predict the critical point and to some extent the elastic 
behavior of the chain in the pinned state. Fig. [2] depicts 
the variation of elastic force (/ e ;) with applied external 
force (/) for different p values. Now, when the string is 
just depinned, the element facing the maximum pinning 
force f™ ax (— 1) must be depinned. So, the condition for 
dcpinning is 

f+fei = r; Lax = i. (ii) 

From Fig. [2] it is seen that for all p values the elastic force 
increases with external force upto the point when their 
sum become unity. At that point the dcpinning transition 
occurs. As we will elaborate later, independent estimate 
of the critical points agree well with this measure. Now 




FIG. 2: Model I: The clastic force |/ ei | is plotted 
against external force for different values of p. It is seen 
that the response is piecewise linear, with a change of 
slope associated with elimination of zero pinning force 
on the string. Top left inset confirms this statement, it 

also shows a increase in density of zero pinning force 
(Co) after dcpinning. The bottom right inset shows the 
variation of the point where change of slope occurs (ft) 
and also the depinning force f c ). The continuous lines 
are predictions for small p values (see text). For the 
simulation data shown, we have used L = 5000 and 
Euclidean metric for measuring the length of the string. 

to understand the elastic behavior, first note that when 
the string is pinned, the points facing the weakest pinning 



must also be pinned. Of course, the weakest pinning for 
a particular pinned configuration of the string need not 
be the weakest pinning possible (i.e., 0). It will actually 
depend upon the external force. The pinning condition 
is 

/™ m + fel = /• (12) 

Before this condition is reached, both terms in the left 
hand side would increase in magnitude until the con- 
dition is satisfied. This can be understood from the 
observation that the sites facing the minimum pinning 
would move forward and eventually face a stronger pin- 
ning, thereby increasing the first term. The overall elon- 
gation of the chain during this process would increase 
the second term. The minimum value of the both terms 
for which the above condition is satisfied will be when 
/™ m = fel = f /• This is also verified numerically. This 
indicate the slope of the p = line to be 0.5, which is 
indeed the case (see Fig. [2]). 

For the cases with finite p, f™ m = until the string 
attains a configuration where there is no longer a zero 
pinning force in its surface. Until that point, however, 
Eq. (fT2"|) implies f e i = f, which is actually the case seen 
in Fig. [2] Once there is no more zero pinning on the 
string, it behaves as p = case and thereby attains a 
slope 0.5 as before. It is also checked numerically that 
the change of slope is indeed associated with the absence 
of zero pinning in the string (see top left inset of Fig. [2j . 
Furthermore, when p is very small, a small value of / is 
required to move all the sites having zero pinning and to 
reach the elimination of zero pinning point. In the first 
approximation one can assume that essentially the zero 
pinned sites are the only ones to move as f e i will be small 
as well. Under this approximation, the change of slope 
would occur at / = p. As one can see (bottom right 
inset of Fig. [2]), this is valid almost upto p = 0.2. In this 
regime, the critical point would be f c = (2 — p)/3, which 
is indeed the case for small p. For large values of p, there 
cannot be a pinned configuration where there is no zero 
pinned site. Therefore, for large values of p, the f e i = f 
relation is always valid, as can also be seen numerically. 

2. Critical behavior and scaling relations 

Here we present the numerical results of the scaling anal- 
ysis mentioned in III CI for Model I. 

In Fig. O the average velocities are plotted for different 
values of the external force on both sides of the estimated 
critical point f c = 0.500±0.001, which matches very well 
with the estimate a/2 argued above (for p > 0.5). In this 
log-log plot, the power-law decay of the velocity at the 
dcpinning point comes as s straight line from which 5 is 
estimated to be 0.60 ± 0.01. 

Then inn Fig. [3] we plot the saturation values of ve- 
locity in the depinned region. Knowing the critical point 
accurately from the above analysis, the so called velocity 
exponent is estimated to be 0.83 ± 0.01. 



FIG. 3: Model I: The saturation values of the average 
velocities (order parameter) are plotted for different 
external force / = 0.5005, 0.501, 0.502, 0.504, 0.506, 

0.508. The log-log plot gives the estimate for 
6 ps 0.83 ± 0.01. System size is L = 5000. The line is 
guide to the eye. The bottom right inset shows the 
variation of order parameter with external force. The 
top right inset shows the log- log plot of W 2 (t) with t at 
f = f c giving 2/3 = 0.70 ± 0.05. Euclidean metric is 
used for measuring the length of the string. 



In top left inset of Fig. [3] the growth of W 2 {t) at the 
critical point is plotted against time. The power law 
fit gives 2/3 = 0.70 ± 0.05. This growth, however, has 
a bound as in the case of other mean-field models pll j. 
From our estimate S + j3 = 0.95, which is in reasonable 
agreement with the scaling relation. 

In Fig. [5] we plot vt~ s against t\f — f c \ v . We already 
know the value of 8. So by only tuning v we get a data 
collapse and the corresponding estimate for v = 1.35 ± 
0.05. From our previous estimates of 9 and 8 we find 
v = 9 1 '8 ps 1.38, which is in reasonable agreement with 
our independent estimate. 



FIG. 4: Model I: The average velocities are plotted for 
different values of external forces above, below and at 
criticality. In the inset t\f — f c \ u is plotted against vt s 
for different values of /. Estimating <5 (ps 0.60 ± 0.01), 
the fitting value of v is found to be 1.35 ± 0.05. The 
system size is L — 5000 and At — 0.01 for all the 
simulations. Euclidean metric is used for measuring the 
length of the string. 



to 

hi(t + i) =hi(t) + -j- ]T [Mt) - h j+1 (t)\] Ci 

3 

+Vi(hi{t)) + f (13) 

with symbols having usual meanings. 

Note that the force due to extension is still felt equally 
by each element of the string. Therefore, the estimate 
of the critical point, as mentioned before, remains un- 
changed (i.e., half of the maximum pinning force). More 
importantly, the values of the exponents (5, 9, v and 
/3) remains same upto the error bar of their estimates 
(8 = 0.60 ± 0.01, 9 = 0.84 ± 0.01, v = 1.35 ± 0.05, 
P = 0.35 ± 0.05; satisfying very closely the relations 
8 + 13 = 1 and v = 9/8). 



B. Model II 



3. Universality 

We have checked the universality of the above results 
by changing the metric of measuring the distance or the 
elongation of the string. Instead of measuring the length 
using the Euclidean metric as mentioned above, we have 
used the rectilinear length or the Manhattan length of the 
string. In general, the rectilinear or Manhattan distance 
between two points (x\, j/i) and (x-i, y-i) is defined as \x\ — 
1 H— 1 2/i — 2/2 1 - Using this metric, Eq. @ will get modified 



1. Phase diagram 

Consider two sites i and j. Regarding the pinning force, 
let i-th site faces a force which is a local minimum (in 
magnitude) and j-th site is the nearest site having pin- 
ning force at a local minimum. Initially, j-th site will 
move forward and «-th site will lag behind. Therefore, 
the forces on i-th and j-th site will respectively be 

F l = fel - S\ + / 

F j = -fel -H + f, (14) 
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FIG. 5: Model I: The average velocities are plotted for 
different values of external forces above, below and at 
criticality for Manhattan distance metric. In the inset 
t\f ~ fc\ v is plotted against vt s for different values of /. 
Estimating 8 (« 0.60 ± 0.01), the fitting value of v is 
found to be 1.35 ± 0.05. The system size is L = 5000 
and At = 0.01 for all the simulations. 



where /" denotes pinning force at a-th site. The above 
equation indicates that it is only possible for the z-th site 
to come in the same level as that of j-th site is when 

2fJ > Af p , (15) 

where Af p denotes that difference between the pinning 
forces at the two sites. Again, the i-th site can move 
forward only when 

f l p < fel + /■ (16) 

The above two conditions together imply 

fel > f- (17) 

As we shall see later, this condition prevents independent 
ballistic motion of the birds; this is also clear from the 
fact that f e i is a kind of nearest neighbor interaction 
that invokes cooperitivity and the width of the profile 
increases sublincarly with time in the depinncd state (we 
call this weakly depinned state; see Fig. [6]). 

Note further that the maximum pinning force a site 
may face is /™ aa: . Therefore, the condition for depinning 
would be 

f p naX < + /, (18) 

with the equality determining the critical points. 

Now consider the case when / > f e i. Irrespective of 
the magnitude of /, this is a depinned state. This is 
because, as the lower limit of the pinning force is zero, 
for any value of /, some sites, which have pinning less 
than /, will be moving. Now, as / e ; < /, those sites will 




✓ 



/' depinned 

(strong) 



FIG. 6: The phase diagram for Model II depicting 
pinned phase and two types of depinned phase. The line 
ABC follows the equation f el +f = f™ ax , the line OB 
is fel = /■ The strongly depinned phase is where 
Eq. (|17[) is satisfied. Here each site moves forward 
independently. Therefore ballistic motion is seen giving 

linear increase in the width of the profile. The 
boundary between pinned and strongly depinned states 
does not represent a phase transition as in the later 
phase no cooperative phenomenon takes place. The 
weakly depinned region is characterised by sublinear 

growth of the width of the profile. The transition 
between pinned and weakly depinned states is a phase 
transition whose exponents are reported (see text). 



never stop. This is, therefore, a trivial case of depinning 
where the particles move essentially independently and 
the 'width' of the profile increases linearly with time; we 
called this strongly depinned state (see Fig. [5]). 



2. Critical behavior and scaling relations 

Here we report the numerical estimates of the exponent 
values for the pinned and weakly depinned transition. In 
Fig. [7J the order parameter (measured in the rest frame 
of the slowest moving site) is plotted against / — f c to 
estimate the order parameter exponent. It is found to 
be 9 = 1.00 ± 0.01. The width of the profile increases in 
a power law with exponent (3 = 0.65 ± 0.01. It is also 
shown that in the strongly depinned region, the width 
increases linearly. Another quantity is measured here, 
which is the concentration of the sites having positive 
(open upwards) curvature (z c ). This quantity shows non- 
trivial power-law relaxation only for this model with an 
exponent 0.278 ±0.001. 

In Fig. [51 the relaxation of the order parameter in time 
is shown near a critical point. At the critical point, the 
velocity decays in a power law with an exponent 8 ~ 
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0.34 ± 0.01. Also, using the scaling techniques discussed 
in III CI the exponent v is estimated to have the value 
2.95 ±0.05. 

As one can see from the above discussions, the scaling 
relations 6 + 6 = 1 and 6 = 6/v arc almost satisfied. 
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FIG. 7: Model II: The top left inset shows the growth of 
the width for strongly (upper curve) and weakly (lower 
curve) depinned states. In the strongly depinned state, 
the width increases linearly and in the weakly depinned 

state it increases sublinearly (with an exponent 
f3 = 0.65 ± 0.01). The main figure shows variation of the 
order parameter in the pinned and weakly depinned 
transition. It gives 6 = 1.00 ± 0.01 as the order 
parameter exponent. The bottom right inset shows the 

decay of the density of sites (z c ) having positive 
curvature at the pinned and weakly depinned transition 

point. This gives a power-law decay with exponent 
0.278 ± 0.005. This feature is unique for Model II. For 
simulations system size was L = 10 4 . The data 
presented here is for the critical point 
fc = 0.4, = 0.6, f 7 p nax = 1. For other points along the 
phase boundary (between pinned and weakly depinned 
state) the exponents are same. 



IV. SUMMARY AND DISCUSSIONS 

In the present work, we study numerically the scal- 
ing properties of the dynamics of elastic string driven 
through a disordered medium. We have considered two 
models here, in one case the elastic force is of infinite 
range and in second it is short range and constant. 

By infinite range we mean that the elastic force in an 
element depends upon the entire elongation of the chain 
(see Eq. Q). Of course, the sign (direction) of the force is 
not uniform for all the elements; it depends upon the sign 
of the local curvature. We have employed both Euclidean 
metric and Manhattan metric for determining the string 
length. Apart from the elastic force, we have assumed 
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FIG. 8: Model II: The time dependence of the order 
parameter (velocity) is plotted for different values of / 
around f c (at the same critical point mentioned in the 

previous figure) measured in the rest frame of the 
slowest moving site. At f c the velocity decays with an 
exponent S = 0.34 ±0.01. The inset shows the data 
collapse to find the exponent v = 2.95 ± 0.05. The 
simulations arc for L = 10 4 and the estimates of the 
exponent values remains same along the pinned- weakly 
depinned phase boundary. 



a random pinning force and an external force which is 
applied to depin the chain. Using the non-locality of 
the elastic force, we have argued that the critical force 
at which the chain will be depinned is f c = (2a — p)/3 
for small values of p, where a is the maximum pinning 
force and f c = a/2 for p > 0.5. By knowing the critical 
force (for p > 0.5), we integrate the equation of motion 
for the height profile (see Eq.([5])) and perform standard 
scaling analysis to get the critical exponents. In partic- 
ular, we have found the power-law decay of the average 
velocity at the critical point (giving 6) , saturation values 
of the velocities in the depinned region gives 9. From 
off-critical scaling the exponent for the diverging tempo- 
ral correlation length (y) is found. The rms width grows 
in a power-law at the critical point giving the exponent 
j3. All these exponent values fit well with the scaling re- 
lations S + (3 = 1 and v = 6/5. The estimates remain 
consistent within error bars when the same is done for 
other values of p as well as by changing the metric of 
the distance measurement. As discussed already, /3, 6, v 
and S values are 0.35 ± 0.05, 0.83 ± 0.01, 1.35 ± 0.05 and 
0.60 ± 0.01 respectively for Euclidean metric (see III C| 
and 0.35 ± 0.05, 0.84 ± 0.01, 1.35 ± 0.05 and 0.6 ± 0.01 
respectively for Manhattan metric (see IIII A 3|) . These 
are summarised in Table 1 and are also compared with 
the other universality classes. Note that the estimated 
exponent value of 6 w 0.83 ± 0.01 in our model is very 
close to the experimental value of (9 f» 0.8±0.1by Ponson 

sin. 
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TABLE I: The exponents values for other known universality classes are compared with the ones presented here. 



Models 


P 


e 
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EW [23] 


0.85 ±0.03 


0.24 ±0.03 


1.73 ±0.04 




KPZ [23] 


0.67 ±0.05 


0.64 ±0.12 


1.35 ±0.04 




l/r z [20] 


0.495 ± 0.005 


0.625 ± 0.005 


1.625 ± 0.005 




MH [11] 


0.841 ± 0.005 


0.289 ± 0.008 


1.81 ±0.1 


0.160 ±0.005 


Model I 


0.35 ±0.05 


0.83 ±0.01 


1.35 ±0.05 


0.60 ±0.01 


Model II 


0.65 ±0.05 


1.00 ±0.01 


2.95 ±0.05 


0.34 ±0.01 



We have also studied a model (as mentioned before) 
where the elastic force is of constant magnitude for each 
element, although its sign (as before) depends upon the 
sign of the local curvature. In this model we have taken 
the pinning force as constant for each element during its 
entire motion. This model is relevant for the study of 
flying flock of birds. As one can see, the pinned force 
may then represent the (un) fitness of the birds and will 
remain more or less constant. The 'elastic force' in this 
case does not depend upon the entire chain of birds, but 
is only of constant magnitude with the sign being de- 
termined by the nearest neighbors (local curvature). As 
one can see, the zero velocity for a flying flock of birds 
has no meaning. But here we consider the velocity of 
the flock from the rest frame of the slowest moving bird. 
The 'depinning' transition in this case implies that in 
the 'depinned' phase, even from the rest frame of the 
slowest moving bird, the 'width' of the flock increases 
with time as opposed to its remaining constant in the 



'pinned' phase. As one can see from Fig. [51 three phases 
are present. In the pinned phase, all birds finally have the 
velocity of the slowest moving one, thereby keeping the 
width of the profile constant. In the strongly depinned 
case, these is no sense of cooperitivity (see Eq. (|T7|) ) and 
the birds move independently, making the width increase 
linearly (see Fig. [7]). In the weakly depinned phase, how- 
ever, the width increases sublinearly. The other expo- 
nents of the this transition (pinned to weakly depinned) 
is summarised in Table 1 and also are compared with 
other universality classes. As one can see, the scaling 
relations 6 + = 1 and 6 = Bjv are satisfied within the 
error bars of the estimates. 

In summary, we have studied two models of depinning 
transitions for 'elastic' strings driven through disordered 
(quenched) media. The exponents estimated do not fall 
into known universality classes. As argued before the 
models (I & II) seem to be relevant for fracture front 
propagation and flocking of birds respectively. 
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